Extended s-wave pairing originating from the aig band in Naa;Co02 ■ y H2O 



in 
o 
o 

(N 

< 
(N 



Kazuhiko Kuroki^, Seiichiro Onari^, Yukio Tanaka^, Ryotaro Arita'^^^, Takumi Nojima^ 
^ Department of Applied Physics and Chemistry, 
The University of Electro-Communications, Chofu, Tokyo 182-8585, Japan 
Department of Applied Physics, Nagoya University, Nagoya, 464-8603, Japan 
^ Max Planck Institute for Solid State Research, Heisenbergstr. 1, Stuttgart 70569, Germany 
Department of Physics, University of Tokyo, Hongo 7-3-1, Tokyo 113-0033, Japan 

(Dated: February 2, 2008) 

Motivated by recent experiments and band calculation results, we investigate the possibility of 
superconductivity originating from the aig band of Naj;Co02 • yH20 assuming the absence of e'g 
pockets. We adopt a single band U-V model that takes into account the local minimum of the aig 
band dispersion at the F point. Using the fluctuation exchange method, we find that an s-wave 
pairing with a gap that changes sign between the inner and the outer Fermi surfaces is enhanced 
due to the coexistence of spin fluctuations at wave vectors bridging the two Fermi surfaces and the 
charge fluctuations near the K point. Possible relevance to the experimental results is discussed. 
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The pairing symmetry and the mechanism of 
superconductivity (SC)lpyj in a cobaltate Na:rCo02 • yH20 
has attracted much attention lately. There have been 
debates mainly on three issues concerning this mate- 
rial: whether the relevant band is aig or e'g, whether 
the spin fluctuations are ferromagnetic or antiferromag- 
netic, and whether the pairing occurs in the spin-singlet 
or in the triplet channel. In the previous studies, we have 
proposed a possibility of ferromagnetic-fluctuation medi- 
ated triplet /-wave pairing originating from the e'g pocket 
Fermi surfaces (FS),0, U whose presence was first pre- 
dicted in the band calculation by Singh, and has been 
an issue of considerable debate in the band calculation 
studies that have followed. H 0, 0, 0, IS E3 There have 
also been some other theories where the e '„ b ands play 
an important role in the occurrence of SC. Il2t 

However, some recent experimental results seem to 
contradict with the e^-triplet-pairing scenario. For in- 
stance, in the angle resolved photoemission (ARPES) ex- 
periments, the e'g bands lie below the Fermi level. 0, 
llEEEE^ Also, l/TiT and Knight shift measurements 
suggest that the spin fluctuation is not (or at least not 
purely) ferromagnetic in the SC samples. [isl Eol E] | 
and a recent theoretical calculation shows that this is 
consistent with the absence of e'g pockets. 22] As for the 
singlet-triplet debate, recent experiments find a decrease 
of Knight shift belo w Tr for magnetic fields applied paral- 
lel to the planes. |2fll2lj Also, the effect of the impurities 
on Tc is found to be small, suggesting a more or less 
isotropic pairing. |2^ |24] Thus, there is now an increas- 
ing possibility of sin glet pairing, although triplet pairing 
is still not ruled out. [23 

Under these circumstances, here we investigate the 
possibility of SC within the aig band assuming the ab- 
sence of e'g pockets. There already exist some theories 
that focus on the aig band [H Hi IH IH US HI HI 
but our study differs from those at least in one of the 
following senses, (i) We take into account the local mini- 
mum structure of the aig band dispersion at the F point 



found in first principles calculations for bilayer-hydrated 
(BLH) systems, (ii) We adopt band fiUings cor- 

responding to the Co valence (Vco) of ~ +3.5 because 
experiments show that the actual Vco is much less than 
the nominal value ~ 3.65 due to the pres ence of HaO"^ 
and/or oxygen deficiencies. |3J, |33, |3a, |33, 123 Combined 
with (i), an inner FS appears in addition to the outer FS, 
which results in a situation where "disconnected FS" [s^ 
exist. (Note that the Fermi level in the existing band 
calculations for BLH systems corresponds to Vco ~ 3.67, 
so that the inner FS does not exist.) (iii) The nearest 
neighbor (NN) repulsion V, which can be important for 
systems away from half- filling, is taken into account. 

We apply the fluctuation exchange (FLEX) method 
|4oll4ll | to a single band U-V model to obtain the Green's 
function, and then solve the linearized Eliashberg equa- 
tion. We find that the inner FS near the F point strongly 
affects the spin fiuctuations, where the spin susceptibility 
takes its maximum at a wave vector that bridges the in- 
ner and the outer FS. Such spin fluctuations enhance SC 
in a manner that the gap changes sign between the two 
FS. The simplest among such pairing is the extended s- 
wave pairing whose gap does not change sign within each 
FS. This extended s-wave pairing is further enhanced by 
charge fluctuations near the K point induced by the NN 
repulsions. We conclude that a combination of the dis- 
connected FS and the coexistence of spin and charge fluc- 
tuations may lead to a Tc of order O.OOlti (ti is the NN 
hopping) for the extended s-wave pairing, which roughly 
corresponds to the actual Tc of the cobaltate. 

The model is given as H = J2i.jJ2a^v4a(^3o- + 



where we consider 



hopping integrals up to third nearest neighbors, i.e., ii, 
^2, and ^3. ^1 = 1 is taken as the unit of the energy and 
the on-site repulsion is fixed at [/ = 6 throughout the 
study. We mainly focus on the band filling(=number of 
electrons/number of sites) of n = 1.5, which corresponds 
to Vco = 3.5. In the present study, we mainly consider 
the case of ^2 = —0.35 and ^3 — —0.25 (although some 
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FIG. 1: Left panel: The energy dispersion for the non- 
interacting case with t2 — —0.35 and — —0.25. The hor- 
izontal line shows Ef for n = 1.5. Right panel : FS for 
n = 1.5. 



other values are also adopted for comparison), where the 
non-interacting {U = V = 0) band dispersion has a local 
minimum at the F point with a "depth" of 0(10%) of the 
total band width 1421 as seen in FigUI so that an inner 
FS is present for n = 1.5. 

FLEX has been app lied to the Hubbard model with 
on-site U, 0j Q| and has further been a ppli ed to 
the extended Hubbard model by Esirgen et aL[4j. We 
first obtain the renormalized Green's function, G, tak- 
ing bubble and ladder diagrams as the self energy. Here, 
we define the "Fermi surface" by ej^ + Re[E^] — ^ = 0, 
where £j^.is the non-interacting band dispersion, the 
self-energy at the lowest Matsubara frequency and fi 
the chemical potential. We then calculate the pairing 
interaction mediated by spin and charge fluctuations, 
and plug that into the linearized Eliashberg equation, 

AA(fc) = J2k' r(fc-fc')G(fc')G(-fc')A(A:'), where A is 
the gap function, G the Green's function, and F the pair- 
ing interaction with k = (fc,a;„). is the temperature 
where the largest eigenvalue A reaches unity. When the 
NN repulsion V is introduced, the Green's functions and 
the susceptibilities become 7x7 matrices for the trian- 
gular lattice. In the actual calculation, we take up to 
N = 64x64 sites and -(2iVc~ l)7rT < uj„ < {2Nc-1)ttT 
with Nc = 2048. 

As for the pairing symmetries, we focus on s-wave {Ai 
symmetry), d-wave {E2), and /-wave {B2) since these are 
found to be the only pairings that have a chance to dom- 
inate in our model. When V — Q in particular, /-wave 
also has no chance of dominating, since antiferromagnetic 
spin fluctuations alone cannot mediate triplet pairing, so 
we consider only d- and s- waves. 

Let us move on to the calculation results. We first show 
results without the NN repulsion V. When t2 = H = 0, 
the spin susceptibility broadly peaks around the wave 
vectors that bridge the opposite sides of the FS as seen 
in Fig|2Ia). Both d-wave and s-wave pairings have very 
small values of A (Figl^), suggesting that there is no SC 
instability for n — 1.5 at least in the realistic temperature 
range. When we turn on the distant hoppings t2 = —0.35 
and t'i = —0.25 so that the inner FS appears, the spin 
susceptibility still has peaks at the Brillouin zone edge 
^Qs2)i but takes larger values at wave vectors Q^^ that 




FIG. 2: (color online) Contour plots of the spin susceptibility 
for (a) t2 = t3 = and (b) t2 = -0.35 and = -0.25. (c) 
Contour plot the s-wave gap function for = —0.35 and = 
—0.25, where the red solid curves are the FS and the yellow 
dashed lines and curves are the nodes of the gap function. 
(7 = 6, = 0, and T = 0.01. The gaps and the susceptibilities 
are plotted for the lowest Matsubara frequency throughout 
the study. 
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FIG. 3: (color online) A plotted as functions of temperature 
for s-(blue lines with squares) and d- waves (green lines with 
circles). The solid (dashed) lines are for t2 = —0.35 and 
t-i = -0.25 (t2 = = 0). [7 = 6 and V = in both cases. 



bridges the outer and the inner FS as seen in Fig[2Ib). A 
repulsive pairing interaction mediated by the spin fluctu- 
ations around Q^^ can enhance SC because a sign change 
of the gap between the inner and the outer FS can take 
place without adding additional nodes on the FS. In fact, 
A of both extended s-wave (gap shown in Fig[5fc)) and 
d-wave (gap not shown) pairings are enhanced by a factor 
of 3 to 4 at low temperature as seen in Fig|31 Neverthe- 
less, these values are still much smaller than unity even 
at T = 0.0025, which roughly corresponds to the actual 
Tc of the cobaltate. 

This is where the NN repulsion V sets in. The FLEX 
results in the presence of V show that while the spin 
susceptibility Xs has structures similar to those in the 
absence of V , the charge susceptibility Xc peaks near 
the K point (Fig^, which bridges the opposite sides 
of the outer FS. This is because the Fourier transform 
of the NN repulsion takes a large negative value at K 
point. Such charge fluctuations work cooperatively with 
the spin fluctuations around Q to enhance the extended 
s-wave pairing having the same gap sign within each FS 



3 




FIG. 4: (color online) Upper panel: contour plots of the 
spin (xs) and the charge (xc) susceptibilities. Lower panel: 
contour plots of the s- (As), d- (Ad), and /-wave (A/) gap 
functions, where the red solid curves are the FS and the yellow 
dashed lines are the nodes of the gap function. U — 6, V = 2, 
t2 = -0.35, U = -0.25 and T = 0.025. 
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FIG. 5: (color online) A plotted as functions of temperature 
for s-(blue lines with squares), d- (green lines with circles), 
and f-wa.ve (red lines with triangles) pairings. The upper 
(lower) lines are for = 3 (V = 2). U = 6, ^ -0.35, 
and ts = —0.25. A for the /-wave {Bi symmetry) pairing 
in a model for the e'g band is shown for comparison (dashed 
curve). Inset: n dependence of the s-wave A with V = 2 and 
T = 0.025. 



(As in Fig^. This is because the charge fluctuations 
mediate attractive pairing interactions as opposed to the 
spin fluctuations. 

In fact, the eigenvalue of the s-wave pairing is strongly 
enhanced in the presence of V as seen in Fig|Sl d- and 
/-wave pairings are also enhanced because these gaps 
can exploit the attractive interaction near the K point 
at some portions of the FS,^ but s-wave pairing dom- 
inates over d and / at least for the parameter values 
adopted here. Although we cannot go down to low tem- 
peratures in the presence of V, the nearly linear increase 
of log(A) with decreasing log(r) suggests that A may 
reach unity at least around T — O(O.OOlii) for the s- 




FIG. 6: (color online) A plot similar to Fig|^with t2 — —0.15 
and ts — —0.25, V = 2. inset: the energy dispersion for 
U = V = 0, t2 = -0.15, ts = -0.25, [/ = 6, V = 2, and 
T = 0.025. 



wave pairing. Note that the values of A for the s-wave 
pairing is even larger than that for the /-wave pairing 
(Bi symmetry) state originating from the e'g pockets an- 
alyzed in a single band model in ref.'3|, at least for the 
temperature range where the calculation is performed. 

The inner FS in FigQl is smaller than that in Fig^ 
which means that the correlation effect makes the depth 
of the local minimum of the band shallower, so that the 
Fermi level is located very close to the local bottom. We 
find that the values of A are not sensitive to the choice of 
the values of the distant hoppings as far as this relation 
between the Fermi level and the energy at the F point 
does not change significantly. For instance, if we take 
t2 = —0.45, t3 — —0.1, and — —0.04, where t4 is the 
fourth nearest neighbor hopping, the relation between 
the Fermi level and the local minimum at F is similar to 
that for t2 = -0.35 and ts = -0.25, but the outer FS 
has a more hexagonal shape rather than the round one 
shown in FigQ] Performing FLEX calculation for this set 
of parameter values, we find that the values of A are very 
close to those shown in Fig|Sl 

On the other hand, the situation changes drastically if 
the local minimum is too shallow to have an inner FS at 
n = 1.5, as expected from the results for = (FigO. 
We show in Fig|Hlthe calculation results for t2 = —0.15 
and t3 = —0.25, where the local minimum of the band at 
the F point lies far above the Fermi level for n ~ 1.5 as 
seen in the inset. It can be seen that the s-wave state is 
strongly suppressed compared to the case of t2 — —0.35, 
and none of the pairing symmetries are likely to have a 
Tc in the realistic temperature range. 

We further find that when the local minimum is too 
deep (by taking, for instance, a too large — ^2), so that 
the inner FS becomes large, the s-wave pairing is again 
suppressed. This is because the distance between the two 
FS becomes short, so that the gap nodes are located too 
close to the FS. Such a tendency can also be captured 
by fixing t2 = —0.35, is = —0.25, and varying n, where 
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we find that A is maximized around n = 1.5 (the inset 
of FiglSJ. These results show that SC can be reaHzed 
(in the reahstic temperature range) only when the Fermi 
level is located near the local bottom of the band. 

Finally, let us discuss the relevance of the present 
results to the experiments, (i) SC, accompanied by 
enhanced non- ferroma gne tic (or incommensurate) spin 
fluctuations [l^ Ell l45l|. is observed only in the BLH 
samples. 5^ In the present view, this can be explained as 
follows. In the non-hydrated and MLH systems, there is 
a large three dimensionality in the aig bands, as can be 
seen from the split of the two bands (originating from two 
C0O2 layers in a unit cell) at the F point. |E0E3| This 
results in one band without the inner FS and the other 
with 3D inner and 2D outer FS, which would smear out 
the Qsi spin fluctuations discussed here. Moreover, for 
non-hydrated systems inparticular, the local minimum 
depth tends to be small, [3 so that there may be no inner 
FS at all. (ii) s-wave pairing is consistent with the ab- 
sence of time reversal symmetry breaking suggested from 
^SR experiments. ^3>S| (iii) s-wave pairing is also con- 
sistent with the decrease of Knight shift below T.~.|2oLEH 
(iv) Although the nodes of the gap do not intersect the 
FS in the extended s-wave state, the change of both the 
sign and the absolute value of the gap between the inner 
and the outer FS may explain some of the unconventional 
behavior observed in 1/Ti 0,113,01 or heat capacity 
measurements. [2^. ISOj In fact, we have preliminary re- 
sults showing that the absence of a coherence peak and 
the crossover from ~ T^^ to ~ T behavior in 1/Ti can be 
understood within the present s-wave pairing by taking 
into account the effect of impurities, (v) The extended 
s-wave pairing r nay also account for the small impurity 
effect on rc|23j uM since ^j^^pg A(fc) = is not sat- 
isfied due to the asymmetry between the two FS. (vi) 
It has rec ently been found that Tc is maximized around 
^Co ~ 3.5.(33 This is consistent with our conclusion that 
SC can be realized only when the Fermi level is located 
near the local bottom of the band (inset of FigElnote 
that Vco = 5 — n). 

To summarize, we have shown that an extended s-wave 
pairing can arise from the aig band due to a combination 
of disconnected inner and outer FS and coexisting spin 
and charge fluctuations. The present scenario seems to 
be consistent with many of the experimental results. 
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